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SUMMARY 

Near-surface  measurements  of  infrared  (8  to  12  jfu^horizon  radiances  are  pre¬ 
sented  for  2  days  with  differing  surface  wind  speed  conditions  and  vertical  profiles  of 
meteorological  parameters.  A  representative  aerosol  size  distribution  model  for  each 
day  was  chosen  wl^ch  allow^  calculationa^jising  LOWTRAN  6  (Kneizys,  et  al.,  1983) 
match  the  measured  horizon  pixel  radiance  and  the  atmospheric  optical 
depths  d^tehnined  from  the  upwelling  solar  radiances  detected  by  Channel  1  (0.58  to 
0.65(^^)of  the  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  on  board  the 
Natiraal  Oceanic  and  Atmospheric  Administration  (NOAA)-9<^salellite. 

The  selected  models  were  then  used/iercaleulations^'we  apparent  sea  surface 
radiance  as  a  function  of  zenith  angle  using  a  modified  version  of  LOWTRAN  6.(^ol- 
lenweber,  1988).  This  program  incorporates  a  statistical  wave  slope  model-(Cox  and 
Munk,  1954)^  account  lor  the  contributions  to  the  apparent  sea  radiance  from  sly 
reflections  and  emissions  from  the  surface  wave  facets  and  atmospheric  path.  Good 
agreement  between  calculated  sea  radiances  and  those  measured  at  an  altitude  of 
33  m  was  obtained  within  1^^^  the  horizon  for  low  and  moderate  surface  wind 
speeds.  ^ 


■-"4  c. 


-  Using  these  models,  the  relative  contributions  of  the  three  components  to 
radiances  received  at  different  altitudes  from  other  zenith  an^es  were  determined. 
The  calculations  showed  the  emission  of  the  intervening  atmosphere  to  be  the  mcgor 
contributor  to  the  apparent  sea  radiance  received  at  higher  elevations  from  zenith 
angles  within  a  few  degrees  of  the  horizon.  The  resulting  variations  of  apparent  tem¬ 
perature  of  the  sea  stuface  with  zenith  angle  were  also  found  to  be  stron^y  depend¬ 
ent  on  the  vertical  profiles  of  air  temperature.  F  Ihrfrr  o.4ni c-t.  p i  r.i , 
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INTRODUCTION 


To  be  of  practical  use  in  tactical  situations,  airborne  Forward  Looking  Infra¬ 
red  (FLIR)  ^tems  must  be  able  to  detect  a  surface  target  at  ranges  of  several  tens  of 
kilometers.  This  requirement  limits  the  depression  viewing  angles  to  within  a  few 
degrees  of  the  horizontal.  As  an  example,  the  range  to  a  target  as  viewed  with  a  1° 
depression  an^e  is  80  km  for  an  airborne  FLIR  operating  at  an  altitude  of  1  km. 
Depressing  the  viewing  angle  another  2°  places  the  target  only  20  km  away.  To  pre¬ 
dict  the  FLIR’s  performance,  the  spectral  radiance  contrast  between  the  target  and 
the  sea  background  is  required.  The  performance  range  is  determined  to  be  that 
range  where  the  difference  between  the  apparent  ship’s  temperature  (i.e.,  the  actual 
ship  temperature  degraded  hy  the  atmospheric  transmittance)  and  the  effective  back¬ 
ground  temperature  of  the  sea  is  equal  to  the  minimum  detectable  temperature  differ¬ 
ence  (detection)  or  the  minimum  resolvable  temperature  difference  (identification  and 
classification)  of  the  FLIR  system.  For  the  larger  viewing  angles,  the  infrared  (IR) 
radiance  of  the  sea  surface  is  the  primary  contributor  to  the  background  scene.  Closer 
to  the  horizon,  however,  the  sky  radiance  reflections  and  emissions  by  the  interven¬ 
ing  atmosphere  must  be  taken  into  account.  While  a  smooth  sea  is  a  poor  emitter  of 
IR  at  the  small  viewing  angles,  emissions  from  the  individual  wave  facets  of  a  wind- 
ruffled  sea  contribute  to  the  background  scene. 

In  this  report,  the  relative  contributions  of  sl^r  radiance  reflections  and  the 
emissions  from  the  sea  surface  and  the  intervening  atmosphere  are  examined  as  a 
function  of  sensor  altitude.  For  these  calculations,  airborne  measurements  of  the  ver¬ 
tical  profiles  of  meteorological  parameters  and  sea  surface  temperatimes,  obtained  on 
2  days  during  low  and  moderate  wind  speed  conditions,  are  used  as  inputs  to  a  modi¬ 
fied  version  of  LOWTRAN  6  (Wollenweber,  1988).  In  addition  to  the  path  emissions, 
this  version  incorporates  the  Gaussian  distributed  wave  slope  model  of  Cox  and 
Munk  (1954)  to  calculate  the  contributions  from  s^  radiance  reflections  and  the 
emissions  from  a  wind-ruffled  sea  surface.  In  the  following  sections,  the  mathematical 
formulation  of  the  background  radiance  scenes  are  presented.  A  technique  is  also  dis¬ 
cussed  by  which  low-level  measurements  of  horizon  IR  radiances  were  used  with  vis¬ 
ible  upwelling  solar  radiances  (as  detected  by  the  National  Oceanic  and  Atmospheric 
Administration  (NOAA)-9  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  to 
select  the  appropriate  LOWTRAN  6  aerosol  size  distribution  parameters  for  the  radi¬ 
ance  calculations.  The  measured  horizon  radiances  were  also  used  to  test  the  atmos¬ 
pheric  models  used  in  the  radiance  calculations  with  altitude. 


MATHEMATICAL  FORMULATION  OF  BACKGROUNDS 

Consider  the  atmosphere  to  be  composed  of  a  number,  N,  of  isothermal  layers 
characterized  by  temperature  Ti  and  transmittance  T(v,i,(i)  along  the  optical  path 
traversing  the  ith  layer  at  ang^e  )i,  and  v  is  the  spectral  wave  number.  fVom 
KirchofTs  law,  the  radiance  of  the  ith  layer  is 

N{v,i,n)sk  -  [1  -ra(v,/./<)]H'(7,)/ff  (1) 


1 


where  Ta(v,i,)i)  is  the  absorption  transmittance  and  W(Ti )  is  Planck’s  blackbody 
radiation  formula.  Then  the  spectral  radiance  reaching  the  sea  surface  throu^  the 
intervening  atmosphere  is 


/■-I  /-I  (2) 

N(v.i,/t)skl  n  r(v.y,/i)]  =  [l-ra(v./,;<)][  n  T(y,j,n)]W{Tt)/:t. 


Summing  the  contribution  from  all  layers,  the  spectred  radiance  at  the  sea  sur¬ 
face  is  then 


n  i-l  (3) 

N(,V,/i)sk=[  2  [l-ra(v,/.At)][./.|i)](  n  T(v,j./i)]W(Ti)/3t. 

1=1  y=i 

As  shown  in  figure  1,  the  radiance  is  allowed  to  strike  a  wave  facet  on  the 
ocean  surface  with  a  Gaussian  distribution  (Cox  and  Munk,  1954)  of  angular  tilts  a 
and  p  in  the  upwind  and  crosswind  directions,  respectively,  such  that  an  amount 
Niv,ii)'  is  reflected  into  the  sensor  at  an  altitude  Hi  within  the  mth  layer.  The  prob- 
abilily  that  radiance  hits  the  facet  is  equal  to  the  probability  that  the  wave  slope 
exist,  i.e., 

N(v,/iy=P(Sx,Sy)N(v,n)sk  (4) 


where 


PiSx.Sy)  =  l/{2ax(fy)EXP[0.05(Si/<yx*Sj/4)]  (5) 

andSx  =  tantt,Sy  =  tanp,  oj  =  0.003  1.92xl03Vc,  =  3.16x10®  Vc,  with  Vc 

being  equal  to  the  current  wind  speed  in  the  azimuthal  direction  ^  with  respect  to 
the  sensor.  Then,  the  total  spectral  radiance  reflected  from  all  the  wave  facets  into 
the  line  of  sight  of  the  detector  located  in  the  mth  layer  is 

m  - 1 

^(v,®)rj*  =  [  n  T(j,e)]Ui(v,a)PiSx.SyN(y,fj.)rsk)  (6) 

7=1  P 

where  /2(v,n)  is  the  complex  reflectivity  of  sea  water  at  the  reflection  ang^e  A.  In  the 
above  equations,  both  ^  and  A  are  implicit  flinctions  of  Sx  nd  given  by  (WoUen- 
weber,  1988) 


cos n  =  (2Sx/A)  cos  O'  cos  (p  +  {2Sy/A)  cos  0'  cos  ^  -  (B/A)  sin  O'  (7) 

cos  A  =  (2Sx/A)  cos  cos  ^  +  ^Sy/A)  cos  O'  cos  ^  +  ( 1  /A)  sin  P'  (8) 


where  A  =  +  5^  +  1,  8  =  5*  +  5|-1  and  0'  is  the  sensor’s  zenith  angle  at 

the  sea-surfoce  reflection  point. 


2 


N(i'.  9). 


FIgurs  1.  Reflection  geometry  from  a  wind-ruffled  sea  surface. 


Similarly,  the  spectral  radiances  emitted  by  the  sea-surface  wave  facets  ) 
and  the  path  radiance  (Np),  which  reach  the  sensor  at  the  zenith  angle  6,  are  given 
by 

m~l 

N{v,e)„  =  [  n  r(v.y.<?)i5:p(5x.s,)[i -/?(v.n)]M'(r„)/»  (9) 

y  =  1  a 

and 

Niv.  =  [  2  [1  -  Ta(v.  /.  fl)]  (  n  e)]WiTi) /7f  (10) 

^=1  ;=i+l 

where  T„  is  the  sea-surface  temperature,  and  again,  the  an^e  ^  is  implicit  in  the 
reflection  ang^e  A. 

Then,  the  total  spectral  radiance  reaching  the  detector  is  the  sum  of  the  three 
components 


Niv,  ff)tot  =  Niv,  e)„k  +  Niv,  e)ss + Niv,  e)p.  dD 

The  total  spectral  radiance  must  then  be  averaged  over  the  response  of  the  FLIR  sys¬ 
tem,  which  in  this  case  is  taken  to  be  the  8-  to  12-^m  wavelength  band.  Subroutines 
have  been  introduced  into  LOWTRAN  6  (Wollenweber,  1988)  to  calculate  the  total 
band  averaged  radiance.  The  reflection  and  zenith  an^es  are  calculated  with  equa¬ 
tions  5,  7,  and  8  corresponding  to  the  incremented  values  of  wave  slopes  in  the  inter¬ 
vals  -  iox.y  <  Sx.y  <  “iox.y.  To  limit  the  number  of  calculations,  the  zenith  angles  are 
divided  into  a  maximum  of  30  classes  (with  the  criterion  that  each  class  should  con¬ 
tain  at  least  10  percent  of  the  probability),  and  the  averaged  angles  for  each  class  are 
then  used  in  the  radiance  calculations. 


3 


MEASUREMENTS  AND  AEROSOL  MODEL  SELECTION 


For  this  study,  a  Piper  Navfgo  aircraft,  equipped  with  Rosemount  temperature 
and  pressure  probes  and  an  EG&G  dewpoint  sensor,  made  vertical  spirals  over  the 
ocean  to  obtain  the  profile  of  temperature,  relative  humidity,  and  pressure  which  are 
required  inputs  to  the  LOWTRAN  6  computer  code  for  calculating  the  sea  and  sky 
radiances.  A  Barnes  PRT-5  radiation  thermometer  was  also  on  board  the  aircraft  for 
measiirement  of  the  sea-surface  temperatures  from  low-level,  constant-altitude 
flights.  The  vertical  profiles  of  temperature  and  relative  humidity,  measured  at  1440 
PST  on  29  September  and  25  November  1988  off  the  coast  of  San  Diego,  California, 
are  shown  in  figures  2  and  3  respectively.  The  measured  sea-surface  temperatures 
were  15.2**C  on  29  September  and  17.6°C  on  25  November.  At  the  time  the  meteoro¬ 
logical  parameters  were  obtained,  measurements  of  IR  (8  to  12  fim)  horizon  radiances 
were  also  made  with  a  calibrated  thermal  imaging  system  (AGA  THERMOVISION, 
model  780)  using  a  2.95*’  field-of-view  lens  with  an  instantaneous  field  of  view  of  0.87 
mr.  The  response  of  the  ^tem  was  determined  by  placing  a  blackbody  of  known  tem¬ 
perature  (±0.1  **0  for  temperatures  <50‘’C)  in  front  of  the  lens  aperture.  The  digit¬ 
ized  video  signal  transfer  function  of  the  ^stem  then  allowed  the  blackbody 
temperatiu-e  to  be  reproduced  to  within  -irO  2‘’C.  For  these  measurements,  the  scan¬ 
ner  was  located  at  an  elevation  of  33  m  on  the  Point  Loma  peninsula  in  San  Diego 
and  was  directed  due  west  over  the  ocean  such  that  approximately  half  of  the  field  of 
view  was  above  and  half  below  the  horizon.  The  measured  radiance  scenes  are  shown 
in  figure  4.  The  temperatures  of  the  different  colors  in  each  scene  are  also  identified 
by  the  color  bars  displayed  on  the  left  which  correspond  to  the  midpoints  of  the  tem¬ 
peratures  printed  above  and  below  each  bar.  The  data  processing  software  of  the  AGA 
system  also  allows  the  vertical  profile  of  the  effective  blackbody  temperature  in  the 
scene  to  be  displayed  on  the  right  side  of  the  thermogram.  In  each  scene,  the  horizon¬ 
tal  cursor  is  situated  on  the  pixel  corresponding  to  the  maximum  temperature 
(20.4*’C  or  3.46  mW/cm^  sr  and  16.5*’C  or  3.23  mW/cm^  sr)  which  is  t^en  to  coincide 
with  the  infrared  horizon.  The  increase  is  temperature  with  altitude  (decreasing  ze¬ 
nith  angle)  in  the  29  September  scene  reflects  the  increase  in  air  temperature  with 
altitude  shown  in  the  figure  2.  Similarly,  the  decrease  in  temperature  with  altitude  in 
the  25  November  scene  can  be  related  to  the  decrease  in  air  temperature  with  alti¬ 
tude.  The  larger  sky  radiance  values  on  29  September  may  also  partially  result  from 
the  lower  relative  humidity  above  200  m  (figure  3),  which  decreased  the  number  of 
activated  aerosols  available  to  scatter  the  radiation. 

These  measurements  can  be  modeled  with  LOWTRAN  6  calculations  to  aid  in 
selecting  an  appropriate  aerosol  model  for  radiemce  calculations  on  each  day.  The 
LOWTRAN  6  aerosol  model  chosen  for  the  calculations  is  the  Navy  Maritime  Aerosol 
Model.  This  model  is  the  sum  of  three  lognormal  size  distributions  and,  in  addition  to 
the  svuface  wind  speeds  (current  and  24  h  averaged)  and  relative  humidity,  requires 
the  input  of  an  air  mass  factor  which  identifies  the  origin  of  the  aerosols  as  either 
marine  or  continental  and  is  allowed  to  range  between  integer  values  of  1  for  open 
ocean  to  10  for  coastal  regions.  Also,  when  an  accurate  surface  visibility  is  available 
as  an  input,  the  model  is  adjusted  so  that  the  calculated  visibility  at  a  wavelength  of 
0.55  fun  is  the  same  as  the  observed  value.  The  air  mass  factor  is  defined  in  terms  of 
atmospheric  radon  content  or  an  air  mass  trajectory  analysis  to  determine  the  time 
the  air  mass  has  been  over  land.  As  neither  of  these  techniques  was  available,  an 
alternate  method  was  used  to  select  an  appropriate  air  mass  factor.  Using  the  current 
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FIgur*  2.  Profiles  of  air  temperature  measured  with  altitude  on  29 
September  and  25  November  1987  off  the  coast  of  San  Diego. 


Figure  3.  Profiles  of  relative  humidity  measured  with  altitude  on  29 
September  and  25  November  1987  off  the  coast  of  San  Diego. 
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Figure  4.  Thermogram  of  near  horizon  Infrared  (8  to  12  fim)  radiances  measured  over  the  ocean  by  an 
AQA  THERMOVISION  Model  780  thermal  imaging  system  on  29  September  and  2S  November  1987. 


and  24-h  averaged  wind  speeds  (Vc  =7.5  m/s  and  V  =  4.4  m/s  for  29  September  and 
Vc  =  2.5  m/s  and  V  =  3.9  m/s  for  25  November)  measured  on  shore  at  the  AGA  site 
and  the  vertical  profiles  of  meteorological  pau'ameters,  LOWTRAN  6  calculations 
were  made  to  agree  with  the  maximum  pixel  radiance  in  the  scene  using  nonunique 
combinations  of  air  mass  factors  and  visibilities.  In  these  calculations  the  meteoro¬ 
logical  profiles  were  divided  into  33  layers  as  allowed  by  LOWTRAN  6.  The  lower  lay¬ 
ers  of  the  profiles  are  also  divided  into  sublayers  containing  the  same  amount  of 
absorbing  and  scattering  material  and  the  temperature  as  the  original  layer.  This 
artificial  layering  has  been  found  necessary  (Wollenweber,  1988)  to  remove  the 
anomolous  dip  (Hughes,  1987)  which  occurs  when  aerosols  are  included  in  the 
LOWTRAN  6  radiance  calculations  for  zenith  angles  close  to  90°.  As  the  AGA  scan¬ 
ner  could  not  be  accurately  plumbed,  the  zenith  an^e  of  the  infrared  horizon  was 
taken  be  0.01°  less  than  the  angle  for  which  the  LOWTRAN  calculations  indicated 
the  refi’acted  ray  path  first  hit  the  earth.  In  these  cases,  the  zenith  angle  correspond¬ 
ing  to  maximum  radiance  is  90.17°.  In  figure  5,  the  solid  lines  represent  the  locus  of 
points  which  allow  the  LOWTRAN  calculations  to  match  the  measured  horizon  pixel 
radiances  with  the  different  combinations  of  air  mass  factors  and  visibilities.  Note 
that  these  calculations  were  made  using  a  modified  current  wind  speed  component, 

Ai  =  which  is  different  from  the  value  published  in  LOWTRAN  6.  This 

modification  was  necessary  to  match  previously  published  measurements  of  IR  slgr 
radiances  and  near-surface  aerosol  size  distributions  (Hughes,  1987)  using  the  model. 
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Hgure  5.  Loci  of  points  of  LOWTRAN  6  caicuiations  with  different  com¬ 
binations  of  air  mass  factors  and  vtsibiiities  which  match  measured  val¬ 
ues  of  IR  horizon  radiances  (solid  lines)  and  satellite  detected  visible 
optical  depths  (dashed  lines)  for  29  September  and  25  November  1987. 
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At  the  time  of  the  IR  radiance  measurements,  upwelling  solar  radiances, 
detected  by  Channel  1  (0.58-0.65  )im)  of  the  AVHRR  on  board  the  NOAA-9  satellite 
as  it  passed  close  to  San  Diego,  were  used  to  determine  the  total  atmospheric  optical 
depths  (^o)-  They  were  determined  using  the  SAIC  (Science  Applications  Interna¬ 
tional  0)rp.)  satellite  radiance  computer  code,  which  makes  use  of  a  direct  linear 
relationship  between  upwelling  radiance  and  the  total  atmospheric  optical  depth 
(Griggs,  1975).  Although  the  vertical  structure  of  the  meteorological  parameters 
which  control  aerosol  growth  were  different,  the  total  optical  depths  for  the  2  days 
were  nearly  identical.  In  a  manner  similar  to  the  IR  radiance  calculations,  calcvila- 
tions  of  total  atmospheric  optical  depth  were  made  to  agree  with  the  measurements 
on  both  days.  In  these  cases,  the  optical  depths  were  calculated  assuming  all  the  aero¬ 
sols  to  be  confined  to  the  mixed  boundary  layer  as  determined  by  the  aircraft  flights. 
The  calculations  of  optical  depths  with  different  combinations  of  air  mass  factor  and 
visibility  are  also  shown  in  figure  5  as  the  dashed  lines.  The  intersections  of  the  solid 
and  da^ed  lines  then  determine  the  best  combinations  of  air  mass  factors  and 
visibilities  to  be  used  in  the  background  calculations.  The  inset  in  the  figure  shows 
the  good  agreement  between  the  inferred  visibilities  and  those  calculated  using  aero¬ 
sol  size  distributions  measured  at  the  lowest  levels  of  the  aircraft  flights.  These 
visibilities  seem  reasonable  since  Los  Coronados  Islands,  located  between  25  and 
30  km  off  the  coast  of  San  Diego,  were  not  visible  to  the  naked  eye  on  either  day. 

CALCULATION  OF  BACKGROUND  RADIANCE  SCENES 

Figures  6  and  7  show  the  comparison  of  the  measured  and  calculated  IR  radi¬ 
ances  for  zenith  angles  within  about  1**  below  the  horizon  using  the  air  mass  factors 
and  visibilities  determined  for  the  2  days  in  figure  5.  In  both  cases,  the  mtyor  con¬ 
tributor  to  the  total  radiance  just  below  the  horizon  is  the  path  emission.  While  the 
reflected  sl^  radiance  and  the  surface  emission  are  small,  their  contribution  to  the 
total  radiance  at  this  low  level  (33  m)  of  observation  cannot  be  neglected.  It  is  inter¬ 
esting  to  note  only  a  small  reversal  of  the  relative  magnitudes  of  the  reflected  sl^r 
radiances  and  surface  emissions  between  the  two  sets  of  calculations,  which  demon¬ 
strates  the  small  influence  of  the  wave  slopes  for  the  moderate  wind  speeds  on  29 
September.  Both  the  calculated  and  measured  total  radiances  on  both  days  are  in 
good  agreement,  which  places  confidence  in  their  usefulness  with  radiance  calcula¬ 
tions  at  other  altitudes  and  zenith  angles. 

Using  the  selected  aerosol  models,  the  contributions  of  the  path,  sea,  and 
reflected  sl^  radiances  to  the  total  background  radiance  were  calculated  as  a  function 
altitude  and  zenith  angle.  In  figures  8  and  9,  examples  are  presented  of  calculations 
for  both  days  with  a  sensor  altitude  of  1000  m.  For  zenith  angles  less  than  about  95°, 
the  mcgor  contribution  to  the  background  is  the  path  emission  with  the  reflected  sl^r 
radiance  being  less  than  10  percent  of  the  total  in  both  cases.  In  figure  10,  the  result¬ 
ing  apparent  blackbody  temperatures  of  the  sea  versus  zenith  angles  for  both  days  are 
compared.  Again,  the  higher  apparent  temperature  for  zenith  angles  near  the  horizon 
on  29  September  results  from  the  path  emission  from  the  warmer  elevated  layers.  The 
rapid  falloff  of  path  emission  with  increasing  zenith  angle  (i.e.,  shorter  slant  paths  to 
earth)  is  the  cause  of  the  decrease  in  apparent  temperature  on  this  day.  In  contrast, 
the  temperature  increase  on  25  November  is  the  result  of  the  increase  in  emission 
from  the  warmer  sea  surface  with  increasing  zenith  angles.  The  relative  contributions 
of  the  three  components  will,  of  course,  change  for  other  altitudes.  Apparent  sea 
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Figure  6.  Comparison  of  the  measured  and 
calculated  IR  radiances  for  zenith  angles 
about  1  **  below  the  horizon  using  the  air 
mass  factors  and  visibilities  determined  for 
25  November  1987  in  figure  5. 


ZENITH  ANGLE  (deg) 


Figure  7.  Comparison  of  the  measured  and 
calculated  IR  radiances  for  zenith  angles 
about  1  °  below  the  horizon  using  the  air 
mass  factors  and  visibilities  determined  for 
29  September  1987  in  figure  5. 


temperatures  versus  zenith  angles  calculated  for  sensor  altitudes  of  200  m,  500  m, 
and  1000  m  are  compared  for  each  day  in  figures  11  and  12.  The  megor  differences 
occur  on  29  September  close  to  the  horizon  where,  at  a  zenith  angle  of  92°,  the  appar¬ 
ent  sea  temperatures  at  200  m  and  1000  m  differ  by  as  much  as  6°C.  Smaller 
differences  in  temperature  (w2°C)  occur  near  this  zenith  angle  on  25  November.  As 
the  nadir  zenith  angle  is  approached,  the  apparent  temperatures  at  each  altitude 
closely  approach  the  measured  sea  surface  temperature. 


DISCUSSION 

The  results  of  this  study  have  shown  that  the  apparent  background  tempera¬ 
ture  of  the  sea  surface  varies  with  the  altitude  and  viewing  an^e  of  the  sensor.  It  is 
also  strongly  dependent  upon  the  vertical  profiles  of  temperature  and  relative  humid¬ 
ity  and  the  actual  surface  temperature.  The  primaiy  contributor  to  the  background 
scene  for  viewing  angles  close  to  the  horizon  is  the  emission  of  the  intervening  atmos¬ 
phere.  For  the  wind  conditions  examined  here,  however,  there  was  little  difference  in 
the  surface  rou^mess  effects. 

The  maximum  detectable  range  of  a  FLIR  ^tem  is  usually  defined  as  that 
range  where  the  actual  temperature  difference  between  a  target  and  its  background 
(Tt  -  T^)  is  degraded  by  the  atmospheric  transmittance,  riR),  to  an  apparent  tem¬ 
perature  difference,  aTq,  equal  to  the  minimum  detectable  temperature  difference  of 
the  system.  This  approach  ne^ects  the  changing  background  scene  with  sensor 
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Flgur«  8.  Calculations  of  contributions  of  the  path,  sea.  and  reflected  sky 
radiances  to  the  total  background  radiance  as  a  fuK^tion  of  zenith  angie 
as  viewed  from  an  altitude  of  1000  m  on  25  November  1987. 


Figure  9.  Calculations  of  contributions  of  the  path,  sea,  and  reflected 
sky  radiarKes  to  the  total  backgrourxl  radiarK»  as  a  function  of  zenith 
angle  as  viewed  from  an  altitude  of  1000  m  on  29  September  1987. 
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PIgtir*  10.  Calculations  of  the  apparent  biackbody  temperature  of  the  sea 
/ersus  zenith  angle  as  viewed  from  an  altitude  of  1000  m  on  29  September 
and  25  November  1987. 
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Figure  1 1 .  Calculations  of  the  apparent  biackbody  temperature  of  the  sea 
versus  zenith  angle  as  viewed  from  altitudes  of  200  m,  500  m,  arxJ  1000  m 
on  29  September  1987. 
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ngura  1 2.  Calculations  of  the  apparent  blackbody  temperature  of  the 
sea  versus  zenith  angle  as  viewed  from  altitudes  of  200  m,  500  m, 
and  1000  m  on  25  November  1987. 


position  as  the  LOWTRAN  calculations  are  iterated  to  determine  the  solution  to  the 
appropriate  value  of  t(jR).  The  correct  approach  is  to  use  the  LOWTRAN  code  to 
directly  calculate  the  sum  of  the  target  and  path  radiances  received  by  the  sensor  at  a 
range  R  as  N(R)up  =N(R  =  0)tT(R)  +  N(R)p  • 

N{R)up  is  then  converted  to  an  equivalent  blackbody  temperature  T{R)up  by 
an  iterative  solution  to  Planck’s  blackbody  formula.  Similarly,  an  equivalent 
blackbody  temperature,  Tj ,  of  the  total  background  radiance  (figure  12)  at  the 
specified  range  is  calculated  and  the  resulting  apparent  temperature  difference, 
AT{Ra)  =  T(R)up  +  T(R)b.  determined.  'The  intersection  of  the  curves  of  ^T{R)a  and  the 
^tern’s  minimum  detectable  temperature  difference  (MDTD)  plotted  versus  range 
will  then  determine  the  maximum  detectable  (MDR)  of  the  target.  In  real-time 
situations,  the  amount  of  computer  time  required  to  determine  the  MDR,  using  an 
iteration  technique,  could  be  reduced  considerably  by  search  routines  at  speciHed 
zenith  angles  at  each  altitude. 
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